Adiabatic transformation as a search tool for new topological insulators: distorted ternary 
Li 2 AgSb-class semiconductors and related compounds 
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We demonstrate that first-principles based adiabatic continuation approach is a very powerful and efficient 
tool for constructing topological phase diagrams and locating non-trivial topological insulator materials. Using 
this technique, we predict that the ternary intermetallic series Li 2 M'X where M'=Cu, Ag, Au, or Cd, and 
X=Sb, Bi, or Sn, hosts a number of topological insulators with remarkable functional variants and tunability. We 
also predict that several III-V semimetallic compounds are topologically non-trivial. We construct a topological 
phase diagram in the parameter space of the atomic numbers of atoms in IA2M' X compounds, which places 
a large number of topological materials presented in this work as well as in earlier studies within a single 
unified topological framework. Our results demonstrate the efficacy of adiabatic continuation as a useful tool 
for exploring topologically nontrivial alloying systems and for identifying new topological insulators even when 
the underlying lattice does not possess inversion symmetry, and the approaches based on parity analysis are not 
viable. 



The recent discovery of topological insulators has opened 
up a new research direction in condensed-matter and materials 
science communities with tremendous potential for practical 
applications as well as for conceptual novelties PES The un- 
usual properties of topological insulators become manifest at 
their surfaces through the appearance of spin-polarized metal- 
lic states even though the bulk remains an insulator. These 
surface states have been predicted to host a large number of 
interesting topological phenomena.®^! For instance, the elec- 
trodynamics of a topological insulator entails additional terms 
in Maxwell's equations yielding novel Faraday-Kerr effect, 
magnetoelectric effect, and axion physics. 56 Other exciting 
possibilities with topological surface states involve exciton 
condensation 7 , presence of an image magnetic monopole in- 
duced by an electric charge, and Majorana fermions induced 
by the proximity effect from a superconductorP2^i n this way, 
topological insulators may be thought of as unifying con- 
densed matter physics with particle physics and cosmology. 4 

The practical realization of many of the aforementioned 
opportunities with topological insulators has not been possi- 
ble to date because the choice of available topological ma- 
terials is quite limited. The currently known topological 
insulators mainly belong to the Z 2 class, which inc ludes 
Bi 2 Te3 famil}I^EI, the ternary half-Heusler compounders 
thallium-based chalcogenided^EU, Ge m Bi 2n Te m+ 3 n -clas^2, 
quaternary chalcogenide compounds, and ternary famatinite 
compounds 20 . The search has also been extended to other 
classes of materials where many-body physics and/or crys- 
talline symmetries conspire to produce non-trivial topological 
phases I2E9 Despite this progress however the realization of 
wide classes of emergent topological properties has remained 
a challenge. 

In connection with the materials discovery effort, adiabatic 
continuation is an elegant and powerful tool for obtaining a 
handle on the Z 2 time reversal invariant even when the lattice 
does not possess a center of inversion symmetry. The idea 
is to start with a known topological insulator and exploit the 



fact that if the Hamiltonian of this system can be transformed 
into that of the new system of interest through a continuous 
series of transformations without inducing a band inversion, 
then the new material will also be a topological insulator. In 
what follows, the adiabatic continuation approach will be seen 
to be not only a method for identifying the topological prop- 
erty of a material, but to also provide a pathway for finding 
new classes of topologically interesting materials. We specif- 
ically discuss how gray Sn, which is known to be topolog- 
ically non-trivial, can be mutated to obtain new families of 
topologically non-trivial phases in zinc-blende-like materials 
and their superlattice variations. In particular, we predict that 
distorted-Li 2 AgSb, a lightweight compound, is adiabatically 
connected to Sn, and thus harbors a 3D topological insulator 
state and that its cubic ground state lies near a critical point. 
In contrast, the related Li2CuSb-type compounds are found to 
be topologically trivial. A topological phase diagram in the 
parameter space of the atomic numbers of the constituent el- 
ements in Li 2 AgSb structure is also adduced, which places 
a number of topological materials within a uniform topolog- 
ical framework. The ternary Li 2 M'X series (with a number 
of variants) is thus identified as a new platform for deriving 
a host of topological compounds, alloys and nanoscale het- 
erostructures via the flexibility of their lattice parameters and 
spin-orbit interaction strength. We also predict several III-V 
semimetallic compounds to be topologically nontrivial. 

We start our discussion of band inversion in the electronic 
structure with reference to gray tin 26 . It is useful to under- 
stand first how a non-trivial topological phase emerges in Sn, 
whereas its isostructural counterparts such as silicon (Si) and 
germanium (Ge) are topologically trivial, see Fig. 1(c). As 
shown in Fig. 1(d), the s-like conduction bands (blue dots) 
lie well above the p-like valence bands in Si and Ge, indi- 
cating the absence of band inversion in Si and Ge. As the 
lattice constant increases upon replacing Si or Ge with Sn, the 
conduction band moves down to touch the valence band, re- 
sulting in a band inversion at only the T -point. This pushes 
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FIG. 1 . (Color online) Crystal structure and topological band inver- 
sion, (a) Crystal structure of Li 2 AgSb. Li, Ag, and Sb are denoted 
by green, red, and blue balls, respectively. Sb and Ag form the zinc- 
blende sublattice. (b) The zinc-blende structure of InSb. In and Sb 
are denoted by gold and blue balls, respectively, (c) The diamond 
structure of gray tin. Diagrams in (d) (e), and (f) illustrate band struc- 
tures near the T-point for trivial insulator, non-trivial semimetal, and 
non-trivial insulator, respectively. Blue dots denote the s-like orbitals 
at T. Band inversion occurs in the non-trivial case where the s-like 
orbitals at the T-point fall below the four-fold degenerate j = 3/2 
states. The degeneracy of the j = 3/2 states is lifted by the lattice 
distortion in the non-cubic case. 



the twofold degenerate s— like bands, which were originally 
lying in the conduction band in Si or Ge, below the valence 
bands of the fourfold degenerate p— states of total angular mo- 
mentum j = 3/2, and the Fermi level (Ep) goes through the 
p— states as shown in Fig. 1(e). As the s— states are now oc- 
cupied at T— point, the Z 2 topological invariant picks up an 
extra time-reversal factor of -1 compared to Si or Ge. There- 
fore, Sn is a topologically nontrivial zero gap semiconductor 
or semi-metal in its pristine phase Pa 

If we define band-inversion strength (BIS) as the energy 
difference between the j = 3/2 states and the twofold degen- 
erate 5-states at the T-point, compounds with positive values 
of BIS will be topologically non-trivial, while those with neg- 
ative values will be topologically trivial. The discussion in 
the preceding paragraph suggests a general route for realizing 
a topologically insulating phase where one lifts the 4-fold de- 
generacy of the j = 3/2 states at the T— point by breaking 
the cubic symmetry of the crystal without inducing a nega- 
tive band-inversion strength. We explore this possibility by 
showing how interplay of effects of finite distortion, heavy el- 
ement substitution, and/or formation of superlattice of cubic 
structure, can either remove the band inversion if the band- 
inversion strength is not sufficiently large, or conversely, it 
can enhance the band-inversion strength. 

In carrying out adiabatic transformations, it is important to 
consider the crystal structure first. Here, we begin with the 
ternary compositions Li 2 M' X (M'=Cu, Ag, and Au, X=Sb 
and Bi) whose crystal structure belongs to the space group 
F43m, with the atomic arrangement presented in Fig. 1(a). 
M' and X atoms occupy the Wyckoff 4d and 4a positions, 
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FIG. 2. (Color online) Electronic band structures of topological 
insulating state in Li2AgSb with lattice distortion, (a) Li2AgSb with 
experimental lattice constants, (b) Li2AgSb with 3% expansion of 
lattice constants, (c) Rhombohedral Li2AgSb with hexagonal lattice 
constants a = ao + 3%ao and c = Co, where ao and Co correspond 
to the experimental values, (d) A sketch of band structure near the T- 
point for topologically non-trivial Li2AgSb with a lattice distortion. 
The s-like states are marked with a blue dot. Lattice distortion causes 
a gap to open at Ef, resulting in a topological insulating state in 
which surface bands will span the bulk band gap, resembling the 
dispersion plotted with red lines. 



respectively. Li atoms fill the remaining empty space in the 
Wyckoff 4b and 4c positions. Because M' and X atoms form 
a zincblende type sublattice, these materials resemble InSb if 
Li atoms are removed [Fig. 1(b)]. Note that In and Sb pre- 
cede and follow Sn in the periodic table and InSb resembles 
gray Sn with diamond structure as shown in Fig. 1(c). These 
observations suggest that Li 2 M'X compounds could be can- 
didates for Z 2 topological insulators if their electronic band 
structures resemble the band structure of structurally similar 
Sn. As already noted above, since there is no spatial inversion 
symmetry in Li 2 M y X or in zinc-blend InSb, parity method^ 
of computing Z2 topological number are not viable. 

First-principles band structure calculations were performed 
with the linear augmented-plane-wave (LAPW) method us- 
ing the WIEN2K package 27 within the framework of the den- 
sity functional theory (DFT). The generalized gradient ap- 
proximation (GGA) 28 was used to describe the exchange- 
correlation potential. Spin orbital coupling (SOC) was in- 
cluded as a second variational step using a basis of scalar- 
relativistic eigenf unctions. 

The computed band structure of Li 2 AgSb along high sym- 
metry lines in the Brillouin zone is presented in Fig. 2(a). 
Away from the T— point, the Fermi level is completely gapped 
and thus the topological properties can be determined from 
observations of band structure only near the r -point. Focus- 
ing on the band structure very close to the Fermi level, we find 
that the orbital angular momentum symmetries of these com- 
pounds are identical to those of the low energy spectrum in 
Figs. 1(d) and 1(e). The s— type and p— type states are nearly 
degenerate at the T-point around which the conduction bands 
and one of the valence band have almost linear dispersion, 
indicating that the system is at the critical point of the topo- 
logical phase [Fig. 2(a)]. With a small expansion of lattice, the 
s/p band inversion occurs as shown in Fig. 2(b). Conversely, 
the topological band inversion in Li 2 AgSb can be removed al- 
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FIG. 3. (Color online) Band inversion strength and topological phase 
diagram. Band inversion strength is plotted as a function of the 
atomic numbers of M' and X elements in the zinc-blende sublat- 
tice. The calculation is carried out at a = 12.8 Bohr which is 3% 
greater than the experimental lattice constant of Li 2 AgSb. The black 
line, which defines the topological critical points by the zero value 
of the band inversion strength, separates the trivial and non-trivial 
topological phases. 



FIG. 4. (Color online) Electronic structure of Li 2 M'X series and III- 
V binary compounds. Bulk electronic structures of TIP in (a), TIAs 
in (b), Li 2 AgBi in (c), Li 2 AuBi in (d), Li 2 CdSn in (e), and Li 2 CuSb 
in (f). The size of blue data points is proportional to the probability 
of s-orbital occupation on the anion site. Grey shaded areas in panels 
(a)-(e) highlight band inversion corresponding to the band- structure 
in Fig. 1(e). 



together by uniformly decreasing all lattice constants, demon- 
strating chemical tunability of the system. 

Since Li 2 AgSb is at the critical point of the topological 
phase, it is possible to achieve a 3D topological phase by a fi- 
nite distortion which lifts the degeneracy of the j = 3/2 states 
and simultaneously induces a band inversion. We find that a 
rhombohedral distortion with expansion along the hexagonal 
ab-plane with a = ao + 3%ao and c = Co, where ao, Co de- 
note experimental lattice constants, does the job as shown in 
Fig. 2(c). The lattice expansion is seen to push the s— states 
well below the p— states, and to also open an insulating gap 
by lifting the four-fold degeneracy of j = 3/2 states through 
the loss of cubic symmetry. Since the band inversion occurs 
only at one point in the Brillouin zone, non-cubic Li 2 AgSb 
is a 3D strong topological insulator. We return to this point 
below to show more rigorously that the aforementioned band 
inversion correctly indicates that distorted non-cubic Li 2 AgSb 
is topologically nontrivial. 

Next, we apply the adiabatic continuity principle to a larger 
range of atomic compositions of Li 2 M'X by systematically 
changing the nuclear charge Z of the atoms. Let us assume 
that atoms at the Li, Ag, and Sb positions possess hypotheti- 
cal nuclear charges Zm = 3 — O.bx + 0.5?/, Zm> = 47 + x, 
and Zx = 51 — y, respectively, where x and y are ad- 
justable parameters. This choice guarantees that the system 
remains neutral for all values of x and y. Note that x and y 
need not be integers and that non-integral nuclear charges are 
easily accommodated in first-principles all-electron computa- 
tions in order to change the Hamiltonian continuously, while 
maintaining charge self-consistency throughout the transfor- 
mation process. This mapping can be started with x = 
and y = 0, which corresponds to Li 2 AgSb, and end with 
x = 3 and y = 1, which corresponds to the artificial com- 
pound He 2 SnSn. He 2 SnSn possess inversion symmetry and 
the wavefunction parity analysis can be used to obtain the 
Z 2 topological invariant 26 . Since He is chemically inert, 
He 2 SnSn has similar band structure as Sn and is a topolog- 



ically non-trivial semimetal with Z 2 = 1. 

A topological phase diagram can now be constructed based 
on the band inversion strength as a function of atomic num- 
bers, as shown in Fig. 3 in the region of < x < 3 and 
< y < 1 at the lattice constant a=12.8 Bohr, which is 3% 
longer than the experimental value for Li 2 AgSb. A borderline 
is drawn at the zero of the band inversion strength (thick black 
line). A phase transition occurs when one crosses this border- 
line in that the two regions separated by the borderline are not 
topologically equivalent. Li 2 AgSb with lattice expansion and 
gray tin are seen to lie on the same side of the thick black 
line, and are thus adiabatically connected, implying that they 
are topologically equivalent. Since Z 2 = 1 for gray tin (actu- 
ally, He 2 SnSn with lattice expansion or distortion), Li 2 AgSb 
with lattice expansion [Fig. 2(b)] or distortion [Fig. 2(c)] is 
also topologically nontrivial with Z 2 = 1. The borderline it- 
self marks topologically critical phases. It would move to- 
ward positive y direction if the lattice constants decreased. 
Li 2 AgSb with experimental lattice constants [Fig. 2(a)] lies 
very close to this critical line. Note that the variation of Zs 
can also be considered as a doping effect in the spirit of the 
virtual crystal approximations poEH Therefore, our topologi- 
cal phase diagram in Fig. 3 could be used also to guide search 
of topologically interesting electronic structures through al- 
loying. 

To further establish the usefulness of the adiabatic con- 
tinuation approach in extending the search for topologically 
non-trivial materials, we carefully examine the phase dia- 
gram of Fig. 3 in which we can identify two new candi- 
dates for topologically nontrivial materials, namely, InSb and 
Li 2 CdSn. InSb is a topologically trivial semiconductor in its 
native phase with the experimental gap of 0.26 eV 32 . Here 
we predict that in InSb with adequate lattice expansion, the 
gap can be closed and a band inversion can occur. It is well 
known that the band gap in the III-V zinc-blende structure de- 
creases for heavier constituent atoms and the lattice constants 
are larger. Earlier first-principles calculations have predicted 
that among heavy-atom III-V compounds, TIP and TIAs are 
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stable in the zinc-blende structures 33 . We take the total-energy 
optimized lattice constants 34 and obtain the band structures for 
TIP and TIAs in Fig. 4. Both compounds possess the band in- 
version feature at the T point. Therefore, thallium-based III-V 
compounds are predicted to be topologically nontrivial. In ad- 
dition, the II- VI compound CdTe lies at x = 1 and y = — 1 
in the phase diagram of Fig. 3. While CdTe is topologically 
trivial, the related compound HgTe is a known topologically 
nontrivial semimetal 35 . Li 2 CdSn is also reported to be a cubic 
crystal 3 ^. If we assume I^CdSn has the same crystal struc- 
ture as Li 2 AgSb, the optimized lattice constant is obtained as 
ao = 12.85 Bohr. Our band structures of I^CdSn in Fig. 4 
show that a band inversion occurs at the T point and the con- 
duction bands are below the Ep around the X -point, allowing 
us to predict that Li 2 CdSn is a topologically nontrivial metal. 

Along the preceding lines, we can theoretically predict 
still other topological insulators. If we replace Li 2 AgSb 
with heavier atoms, the band inversion can be achieved at 
the T— point (Fig. 4). We predict thus that Li 2 AgBi and 
Li 2 AuBi are non-trivial topological metals in their pristine 
phases whereas Li 2 CuSb is a trivial band insulator. Another 
novel way of breaking the cubic symmetry of the zinc-blende 
structure to attain the topological insulating state is to double 
the unit cell, which leads to a tetragonal crystal structure. The 
chalcopyrite semiconductors of I-III-VI 2 and II-IV-V 2 group 
compositions which are described by this tetragonal structure, 
are shown to be non-trivial topological insulators in Ref . [37] 
which is consistent with our predictions. Furthermore, we 
have found that quaternary chalcogenides and ternary fama- 
tinites, which are tetragonal with zinc-blende sublattice, are 
candidates for 3D topological insulators. 2 ^ We note that the 
adiabatic continuation method used in the present calculations 
also predicts that half-Heusler compounds are topologically 
non- trivial.^ 

Concerning the practical realization of topologically non- 
trivial phases in Li 2 AgSb class of compounds, as shown in 
Fig. 2(c), the nontrivial topological insulating phase is pre- 
dicted to exist when a uniaxial tensile strain is applied. This 
could be achieved by growing the materials on suitable semi- 



conducting or insulating substrates. Notably, all the materi- 
als discussed in the present study display electronic structures 
similar to those in HgTe/CdTe, which were used to achieve the 
two dimensional topological insulating phase or the quantum 
spin Hall (QSH) phase in a quantum well system. 35 There- 
fore, our work also expands the list of candidate materials for 
realizing QSH phases. 

In conclusion, we have predicted new classes of topo- 
logical insulator materials in III-V semiconductors and Li- 
based compounds Li 2 M 7 X. Small-gap semiconductors such 
as InSb could be turned into nontrivial phases with a lattice 
expansion. With the example of Li 2 AgSb, we show how 
the nontrivial insulating phase could be obtained by a uni- 
axial strain. The exciting opportunities present within the 
ternary Li 2 M'X series (with a number of variants) offer a 
new platform for realizing multi-functional topological de- 
vices for spintronics and fault-tolerant quantum computing 
applications. Undergirding our study is the demonstration that 
a first-principles approach based on arguments of adiabatic 
continuation provides a powerful tool for materials discov- 
ery in search for new topologically interesting materials en- 
compassing centrosymmetric as well as non-centrosymmetric 
crystal structures. The topological phase diagrams of the sort 
we have obtained for the Li-based compounds in this work, 
which encode global materials characteristics extending from 
single element to binary, ternary and possibly quaternary com- 
positions would provide systematic opportunities for mate- 
rials discovery and tunability between trivial and non-trivial 
topological phases, and serve to motivate further research for 
both fundamental and applied purposes. 
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